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Abstract. The majority of High Mass X-ray Binaries (HMXBs) behave as X-ray pulsars, revealing that they contain a magne- 
tised neutron star. Among the four HMXBs not showing pulsations, and that do not show the characteristics of accreting black 
holes, there is the unusual HMXB 4U 2206+54. Here we present contemporaneous high-energy and radio observations of this 
system conducted with INTEGRAL and the VLA in order to unveil its nature. The high-energy spectra show clear indications of 
the presence of an absorption feature at ~32 keV. This is the third high-energy observatory which reveals marginal evidence of 
this feature, giving strong support to the existence of a cyclotron resonance scattering feature, which implies a magnetic field of 
3.6xl0 12 G. On the other hand, the source is not detected at centimetre radio wavelengths with a 3<r upper limit of 0.039 mjy. 
The expected radio emission for an accreting black hole in the low/hard state, inferred from X-ray flux measurements, would 
be at least 60 times greater than the measured upper limit. Both results firmly indicate that, in spite of the absence of pulsations, 
4U 2206+54 hosts a magnetic accreting neutron star, the first one not to be observed as an X-ray pulsar. 
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1. Introduction 

The vast majority of High Mass X-ray Binaries (HMXBs) 
harbour X-ray pulsars (c.f. Bildsten et al. I1997L believed to 
be young neutron stars with relatively strong magnetic fields 
(B ~ 10 12 G). An important fraction of them are wind-fed sys- 
tems, in which the pulsar accretes from the radiative wind of an 
OB supergiant (perhaps, in some cases, a Wolf-Rayet star). 

Among the handful of HMXBs not displaying X-ray pul- 
sations, only three show the typical characteristics of accreting 
black holes (LMC X-l, LMC X-3 and Cyg X-l). In four other 
HMXBs, pulsations have not been discovered in spite of inten- 
sive searches, but there is no strong evidence identifying the 

Send offprint requests to: P. Blay 
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accreting object as a black hole. In principle, there is no rea- 
son to attribute the lack of pulsations in all these systems to 
any particular characteristic and different models have indeed 
been proposed to explain some of them. There have been sug- 
gestions that 2E 0236.6+6101, whose counterpart is the B0 Ve 
star LS I +61°303, and RX J 1826.2- 1450, identified with the 
06.5 V((f)) star LS 5039, may not be accreting binaries af- 
ter all, but X-ray sources powered by rotational energy from 
a young non-accreting neutron star (Maraschi & Treves IT98 11 
Martocchia et al. 2005 1, although the presence of relativistic 
radio jets points towards the accretion scenario (Massi et al. 
120041 Paredes et al. l2000l In the case of 4U 1700-37, opti- 
cally identified with the 06.5 Iaf+ star HD 153919, a compact 
object of unknown nature and mass M x = 2.4+0.3 M accretes 
material from the wind of the massive supergiant (Clark et al. 

i2om 

The fourth HMXB not displaying pulsations is 
4U 2206+54, identified with the 09p star BD +53°2790 
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(Negueruela & Reig I^OTJTl henceforth NR01). The relatively 
high X-ray luminosity of 4U 2206+54, L x ~ 10 35 erg s" 1 
(at an estimated distance of 3 kpc; NR01), combined with its 
spectral shape, makes the presence of a neutron star or a black 
hole in the system almost unavoidable. There are reasons to 
believe that the compact object in this system is a neutron 
star (NR01; Torrejon et al. 120041 Masetti et al. '2004 1, but the 
possibility of a black-hole has not been ruled out completely 
by previous observations. Analysis of the RXTE/ASM X-ray 
lightcurve revealed a 9.6 d periodicity, which is likely to be the 
orbital period of a compact object (Corbet & Peele 2001 Ribd 
et al. I2005i . Moreover, the X-ray lightcurve displays short 
aperiodic variability, with changes in the flux by a factor ~10 
over timescales of minutes (Saraswat & Apparao l 19921 NR01), 
which are typically seen in wind-fed systems, presumably as a 
consequence of stochastic variability in the wind. 

Interestingly, many high-energy sources not showing pulsa- 
tions are microquasars (containing either black holes or neutron 
stars), while pulsating sources do not show significant radio- 
emission (Fender & Hendry 2000 1. 4U 2206+54 shares many 
characteristics with the well-known microquasar LS 5039 
(Paredes et al. 120001 [2002). Both systems contain a non- 
supergiant late O-type star (Clark et al. 2001) and a compact 
object that does not show pulsations (Ribo et al. Qp99 Reig 
et al. 2003 1 orbiting in a relatively close orbit when compared 
to the majority of HMXBs, and both systems show evidences 
of wind-fed accretion (NR01; McSwain et al.l2004) with X-ray 
luminosities in the range 10 34 -10 35 erg s _1 . However, an in- 
spection of the VLA Sky Survey (NVSS, Condon et al. 1998 1 
reveals no radio emission up to a 3<x upper limit of 1 mJy from 
4U 2206+54. The apparent lack of both radio emission and pul- 
sations does not fit within either the typical scenario of a pulsar 
in a HMXB or the microquasar scenario, and a deeper multi- 
wavelength approach is necessary. 

In this work we present new INTEGRAL and VLA ob- 
servations of the source during the periods 2002 December- 
2004 September and 2003 May-June, respectively. The possi- 
ble presence of a cyclotron line, already suggested by the data 
from RXTE and BeppoSAX (Torrejon et al. 120041 Masetti et al. 
2004 ), and the absence of radio emission are discussed. 

2. Observations and data analysis 

2.1. High-energy observations 

INTEGRAL is a joint European mission in flight from 2002 
October, with three on-board high-energy instruments: the 
Imager on Board INTEGRAL Spacecraft (IBIS), coupled with 
the INTEGRAL Soft Gamma-Ray Imager (ISGRI) and the 
Pixellated Imaging Caesium Iodide Telescope (PICsIT), sen- 
sitive to y-rays from 15 keV up to 10 MeV; the SPectrometer 
in INTEGRAL (SPI), optimised for spectroscopy in the 20 keV- 
8 MeV energy range; and the Joint European Monitor-X (JEM- 
X), consisting of twin X-ray monitors, which provides infor- 
mation at lower energies (3-35 keV). An Optical Monitoring 
Camera (OMC) gives source fluxes in the V (550 nm) band 
and complements the 3 high-energy instruments. All 4 instru- 
ments are co-aligned, allowing simultaneous observations in a 



very wide energy range. A detailed description of the mission 
can be found in Winkler et al. ( 2003 ). 

INTEGRAL observed the region around 4U 2206+54 on 
several occasions during its first 22 months of Galactic Plan 
Survey scans (GPSs), i.e., from 2002 December to 2004 
September. We present in Table ^ a summary of all the 
INTEGRAL revolutions during which the source was inside the 
Field Of View (FOV) of ISGRI. In total, the source was ob- 
served by ISGRI for ~337 ks, but it was significantly detected 
only for 27 ks 1 . For revolutions 70, 74, 145 and 189 only an 
upper limit is given because the source had quite a marginal po- 
sition in the FOV of ISGRI and the detection is not significant 
enough. The JEM-X FOV is smaller and thus data were only 
collected during those revolutions when INTEGRAL pointed 
close enough to the source (marked in Table^with the f sym- 
bol) 2 . Only a few pointings in 2003 May and June fulfill this 
requirement. Although SPI has the largest FOV, it cannot ac- 
quire enough information with one exposure due to the detec- 
tor design. To achieve a S/N~10 for a source like 4U 2206+54, 
SPI would need around 300 ks. Nevertheless, using SPIROS 
in TIMING mode (see Skinner & Connell l2003l a light curve 
has been attained in the 20-40 keV energy range. The obtained 
flux values are in good agreement with the ISGRI data, but 
they have larger uncertainties compared to them. Therefore, 
no data from SPI have been used in this analysis. Data re- 
duction has been performed with the standard Offline Analysis 
Software (OSA) version 4.0, available from the INTEGRAL 
Science Data Centre (ISDC) 3 . A detalied description of the 
software can be found in Goldwurm et al. (2003), Diehl et al. 
( 2003 ), Westergaard et al. ( 2003 1 and references therein. 

Archived RXTEfPCA lightcurves of four long observations 
made between 2001 October 12 and 20 have also been used. 
The Proportional Counter Array, PCA, consists of five co- 
aligned Xenon proportional counter units with a total effec- 
tive area of ~6000 cm 2 and a nominal energy range from 
2 keV to over 60 keV (Jahoda et al. 1996). In order to pro- 
duce lightcurves only the top Xenon layer in standard2 mode 
was used. The durations of these observations range from 17.7 
to 29.8 ks and the complete integration time spans ~100 ks. 
A more detailed description of these observations is given in 
Torrejon et al. (2004). 

2.2. Radio observations 

We observed 4U 2206+54 with the NRAO 4 Very Large Array 
(VLA) at 8.4 GHz (3.6 cm wavelength) on two different 

1 A detection is considered for ISGRI when the detection level, 
which is given by the software package and is not directly related to 
the number of or above background, lies above a value of 8. For the 
typical exposure times of GPS pointings (~ 2 ks), the sensitivity limit 
of ISGRI lies around 2.5 count s at 20-40 keV energy range. 

2 A detection is considered for JEM-X when the detection level, 
which is given by the software package and is not directly related to 
the number of cr above background, lies above a value 20. 

3 http : //isdc . unige . ch/ index . cgi ?Sof t+download 

4 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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Table 1. Summary of ISGRI observations of the field around 4U 2206+54 during the INTEGRAL GPS. We list the revolution 
number, the corresponding dates, the intervals of time when the source was in the FOV of ISGRI during these revolutions in 
Modified Julian Days, the total on-source time within these intervals, the amount of time in which the source was detected, the 
mean count rate and its formal error obtained in the 20^40 keV energy range and the detection level value given by the software 
package (where a significant detection is found for a detection level equal or above a value of 8). For those revolutions for wich 
the significance of the detection was between l<x and 3<x (namely, revolutions 70, 74, 145 and 189) we give 3<x upper limits by 
scaling the detected count rates, assuming that a 3cr detection would correspond to a detection level of 8. The f symbol indicates 
those revolutions in which data from JEM-X allowed the detection of the source. 



Rev. 


Date 


MJD 


On-source time 

(ks) 


Detected time 

(ks) 


Mean count rate 
(count s _1 ) 


Detection level 


26 


2002 Dec 30-2003 Jan 1 


52638.43-52640.07 


14 








31 


2003 Jan 14-16 


52653.32-52655.90 


12 








47 


2003 Mar 03-05 


52701.15-52701.27 


8 








51 


2003 Mar 15-17 


52714.85-52714.96 


8 


2.2 


4.1+0.4 


9.0 


54 


2003 Mar 24-26 


52722.85-52722.93 


6 


2.3 


4.1+0.4 


8.0 


55 


2003 Mar 27-29 


52727.64-52727.67 


4 








59 


2003 Apr 08-10 


52737.03-52737.10 


4 








62 


2003 Apr 17-19 


52746.01-52746.04 


10 


7.3 


3.9+0.4 


8.6 


67f 


2003 May 01-04 


52761.26-52762.45 


8 


6.5 


5.9+0.5 


12.6 


70 


2003 May 10-13 


52769.93-52770.10 


12 




< 7.1 


4.1 


74 


2003 May 22-25 


52781.92-52782.12 


14 


— 


< 4.7 


2.7 


79 


2003 Jun 06-09 


52796.88-52797.07 


14 








82 


2003 Jun 15-18 


52805.94-52806.15 


14 


2.1 


3.2+0.3 


9.5 


87t 


2003 Jun 30-Jul 03 


52820.92-52821.08 


12 


6.5 


5.1+0.4 


11.2 


92 


2003 Jul 15-18 


52835.84-52836.00 


12 








142 


2003 Dec 12-14 


52985.44-52985.60 


12 








145 


2003 Dec 21-23 


52994.41-52994.62 


12 




<7.6 


2.9 


153 


2004 Jan 14-16 


53019.25-53019.45 


14 








162 


2004 Feb 10-12 


53045.49-53045.69 


14 








177 


2004 Mar 26-29 


53090.55-53091.48 


16 








181 


2004 Apr 07-10 


53102.88-53103.05 


15 








185 


2004 Apr 19-22 


53114.56-53114.70 


13 








189 


2004 Apr 30-33 


53126.47-53126.63 


15 




< 4.4 


4.8 


193 


2004 May 12-15 


53138.39-53138.42 


4 








202 


2004 Jun 08-11 


53165.48-53165.67 


15 








210 


2004 Jul 02-05 


53189.33-53189.51 


15 








229 


2004 Aug 29-Sep 01 


53246.96-53247.13 


15 








233 


2004 Sep 09-12 


53258.01-53258.80 


8 








234 


2004 Sep 12-15 


53260.99-53261.92 


17 









epochs: 2003 May 12 from 7:05 to 8:00 and from 11:40 to 
12:52 UT (average MJD 52771.4, during INTEGRAL revolu- 
tion 70) with the VLA in its D configuration, and 2003 May 
20 from 15:27 to 17:20 UT (MJD 52779.7, during INTEGRAL 
revolution 73) with the VLA during the reconfiguration from D 
to A. The observations were conducted devoting 10 min scans 
on 4U 2206+54, preceded and followed by 2 min scans of the 
VLA phase calibrator 2250+558. The primary flux density cali- 
brator used was 1331+305 (3C 286). The data were reduced us- 
ing standard procedures within the NRAO aips software pack- 
age. 



3. Results 

3.1. High Energies 
3.1 .1 . Timing 

Analysis of the X-ray lightcurves clearly shows that the source 
is variable on all timescales. However, except from the 9.6 d 
modulation observed in RXTE/ ASM data (see Corbet & Peele 
2001; Ribo et al. 2005| and believed to be the orbital pe- 
riod, no other periodic variability has been detected so far. 
Unfortunately, the INTEGRAL coverage of the source is not 
enough to test the presence of the orbital periodicity. Therefore, 
orbital period analysis is out of the scope of this paper. 

Pulse period analysis gave negative results for both our 
ISGRI and JEM-X datasets. This was expected, as previous 
searches on similar timescales had also failed (see NR01, 
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Corbet & Peele l200Tl Torrejon et al. 120031 and Masetti et al. 
120051 . 

ISGRI data from consecutive pointings were joined to- 
gether when possible and rebinned to 50 s to search for pos- 
sible longer periods. Nothing was found up to periods of ~1 h. 
The 2CM10 ke V lightcurve and power spectrum for a time-span 
of 6950 seconds during revolutions 67 and 87 can be seen in 
Fig-E The difference between the two lightcurves and between 
the corresponding power spectra is apparent. A quasi-periodic 
feature at ~0.002 Hz (~500 s) can be seen in data from rev- 
olution 87, but it is not present at other epochs. The timing 
behaviour of the source seems to be different in every pointing. 

Little attention has been paid so far to intermediate peri- 
ods (of the order of hours), perhaps because intermediate pe- 
riods would be difficult to detect in the RXTE/ ASM data, es- 
pecially when points are filtered and rebinned as one day aver- 
ages to keep statistical significance. We therefore searched the 
RXTE/PCA lightcurves described in Sect. 12. ll for intermediate 
period pulsations to test the possible presence of a slowly rotat- 
ing NS. Unfortunately, the gaps due to the satellite low-Earth 
orbit are rather large in comparison to the periods searched, 
which certainly hampers somewhat the search. We used epoch 
folding and Lomb-Scargle periodogram techniques with nega- 
tive results. We show in Fig.UJthe power spectrum analysis of 
the whole time span. No significant period is detected, particu- 
larly in the interval [0.95-5.5]xl0~ 4 Hz, which corresponds to 
periods from 3 hours to 30 minutes, approximately. 

Searching inside individual observations, we found some 
periodicities, although with low significance. Inside the 
RXTEfPCA 60071-01-03 observation we find a possible period 
of 6900 s while in observation 60071-01-04 we find 8620 s. 
Folding the entire lightcurve on either of these periods results 
in no pulsed signal. We conclude therefore that the analysis of 
the entire lightcurve does not deliver any significant period in 
the range explored. This result leaves only the possibility of a 
period in the range from some hours to 1 d to be explored. In 
order to do this, a suitable long observation with a high Earth 
orbit satellite, like INTEGRAL, would be required. 

We show in Fig.|3]the long-term lightcurves of 4U 2206+54 
in different energy ranges, from the 2-12 keV of the 
RXTE/ ASM data up to 80 keV for the ISGRI data, spanning 
120 d. As can be seen in the 20^-0 keV lightcurve, an increase 
in brightness occurred during revolution 67 (MJD 52761.36). 
The source brightness increased threefold over a timespan of 
the order of half an hour (see also the top left panel in Fig.Q. 

3.1 .2. Spectral Analysis 

From the whole timespan during which ISGRI collected data 
from 4U 2206+54, the source was inside both the Fully Coded 
Field Of View (FCFOV) of ISGRI and the JEM-X FOV only 
during one pointing in revolution 67 and 3 pointings in rev- 
olution 87. Therefore, the available spectrum from revolution 
67 and the mean spectrum from revolution 87 were used for 
the spectral analysis. Systematic errors of 10% for ISGRI and 
5% for JEM-X were added to our data sets in order to perform 



a more realistic spectral analysis 5 . The software package used 
was xspec 11.2 (Arnaud 1996). 

With the aim of comparing with previously published data, 
the comptonisation model of Sunyaev & Titarchuk ( 1980f, im- 
proved by Titarchukf 1994) including relativistic effects, imple- 
mented in xspec as compTT, and a powerlaw model, modified 
to include photon absorption and a high energy cut-off, were 
chosen to fit the data. For the comptonisation model, the emit- 
ting region temperatures derived were 10+3 and 13 ±8 keV for 
the data of revolutions 67 and 87, respectively. The fits were 
acceptable, with corresponding^ 2 ^ of 1.3 for 173 degrees of 
freedom (DOF) in the first observation, and;^ 2 ed of 1.2 for 176 
DOF in the second observation. The powerlaw parameters of 
both observations are listed in Table[2] 

Both models yield a 4-150 keV flux of ~16xl0" 10 
erg s _1 cirT 2 for revolution 67 and ~8xl0~ 10 erg s _1 cm -2 
for revolution 87. Assuming a distance to the source of 3 kpc 
(NR01), its luminosity amounts to ~1.7xl0 36 erg s _1 and 
~0.9xl0 36 erg s -1 , respectively. Around ~50% of the total lu- 
minosity lies in the 4-12 keV energy band, that is ~8.5xl0 35 
erg s _1 for revolution 67 and ~4.5xl0 35 erg s for revolution 
87. We notice that during these observations the source appears 
brighter than in any previous observation. The RXTE/ASM 
lightcurve confirms that the flux was high in the 2-12 keV band 
as well. 

We show in Fig.|5]the INTEGRAL spectra of 4U 2206+54 
for revolutions 67 and 87. Both spectra suggest the presence of 
an absorption feature around ~30 keV, as already noticed by 
Torrejon et al. (2004) and Masetti et al. (120051 in RXTE and 
BeppoSAX data, see Table|2] An absorption feature through the 
cyclabs model (in xspec notation) was added to the powerlaw 
model. In the revolution 67 spectrum, the absorption feature 
was fitted at an energy of 32+5 keV for a fixed line width of 
3 ke V. The same feature is apparent in the spectrum from revo- 
lution 87, where it was fitted at 32+3 keV 6 . Except for the nor- 
malization factors, the fitted parameters to the datasets of both 
revolutions are compatible between them within the errors (see 
Table|2}. 

There are well known calibration problems in the ISGRI 
Response Matrix Function (RMF) that may cast some doubts 
about the reality of the absorption feature reported here. In or- 
der to investigate if this feature is an instrumental effect, we 
normalised the 20-60 keV spectra of 4U 2206+54 and Crab 
to their respective continua modeled by a powerlaw, and then 
divided the normalised 4U 2206+54 spectrum over that of the 
Crab. We have chosen a Crab observation as close in time as 
possible to our data and with similar off-axis angles, to ensure 
that the RMF and off-axis effects are as much similar as possi- 
ble. We show in Fig.|5]the observed spectra (top panel), their 
ratio to the powerlaw model (middle panel), and the ratio be- 
tween the former ratios (bottom panel). The absorption feature 

5 Private communication from the teams of the instruments. 

6 An F-test was applied to the spectral fits of revolution 87, wich 
has a better signal to noise ratio than that of revolution 67. The im- 
provement of the x 2 by the inclusion of a cyclabs component has a 
12% probability of ocurring by chance. One should take into account 
the limitations of this test when applied to lines (Protassov et al. 2002 1 
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Rev. 87 - 6950 s 
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52820.97 



52821.01 
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52821.05 




0.002 0.004 0.006 0.008 0.01 
Frequency (Hz) 



Fig. 1. ISGRI lightcurves of 4U 2206+54 in the energy range 2CMI-0 keV binned at 50 s (left) and associated power spectra (right) 
for a time span equivalent to 3 INTEGRAL science windows (around 6950 s) during revolutions 67 (top panels) and 87 (bottom 
panels). 




Frequency (Hz) 



Frequency (Hz) 



Fig. 2. Left: Power spectrum of the RXTE/ ASM lighcurve, binned at 25 ks intervals. Note the feature at a frequency of 
^1.2xl0~ 6 Hz, which corresponds to the orbital period (see Corbet & Peele 2001). No other period is found at the limit of 
the lightcurve resolution. Right: Power spectrum of the 100 ks RXTE/PCA lightcurve, binned at 96 s intervals. No significant 
period is detected in the frequency interval [0.95-5.5]xl0~ 4 Hz, which corresponds approximately to the time interval from 3 
hours to 30 minutes. The slope of the power spectrum is ~v~ 0,8 , where v is the frequency. 



around 32 keV is still seen. The quality of the data does not 
allow us to state that the detection is statistically significant, 
but the likely presence of a feature at this position had already 
been reported in the analysis of two other independent datasets 
obtained by two different satellites (see Table|5Jl. As it has been 
seen by three different instruments at different times, the ex- 
istence of this absorption feature is strongly suggested. Such 
features in X-ray spectra are generally attributed to Cyclotron 
Resonance Scattering Features (CRSFs) (see, e.g., Coburn et al. 
2002 and references therein). 

Motivated by the possible presence of this CRSF, we have 
summed up images from those INTEGRAL revolutions with 



significant ISGRI detections. The effective exposure time of 
this mosaic amounts to ~27 ks, and we show its extracted spec- 
trum in Fig. [6] We note that spectral shape changes with lu- 
minosity have been reported in Saraswat & Apparao ( 1992 > 
and NR01. Therefore, by summing up data taken on differ- 
ent epochs we might be losing spectral shape information. 
However, our goal is not to study the shape of the continuum, 
but to achieve an improved signal to noise ratio at the CSRF 
position, which is suggested by both the revolution 67 and rev- 
olution 87 spectra to be at ~32 keV. The best fit to the con- 
tinuum of the new spectrum was a comptonisation model of 
soft photons by matter undergoing relativistic bulk-motion, i.e. 
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Fig. 3. Lightcurves of 4U 2206+54 in different energy bands (quoted inside the different panels). The top panel is for the 3-day 
average of the RXTE/ASM data, where the arrows indicate the epochs of the radio observations. The three bottom panels are for 
the INTEGRAL/lSGRl data in intervals of ~2 ks. Error bars represent the mean of the 1-day average data errors for RXTE/ASM 
data and the formal errors for INTEGRAL/lSGRl data. The two triangles in the 20—40 keV lightcurve indicate 3cr upper limits. 



Table 2. Comparison of most recently published powerlaw+high-energy cutoff parameters, adding photon absorption at lower 
energies, for 4U 2206+54. We remit the reader also to tables 2, 3, and 4 of Saraswat & Apparao ( 1992 1 for earlier EXOSAT 
data. All fluxes shown are unabsorbed in the quoted energy ranges. The errors are at a 68% confidence level for the entry of 
Negueruela & Reig (2001 1 and at 90% confidence level for the rest of the data except for the entry of Corbet & Peele (2001 ), 
where confidence levels were not reported. 



Ref. (Mission, year) 


r 


Ecut 


Eton 


2V H x 10 22 


2 


Flux xlO- 10 


Energy ranj 


;e Boyd 






(keV) 


(keV) 


(atom cm -2 ) 


(DOF) 


(erg s _I cirT 2 ) 


(keV) 


(keV) 


Negueruela & Reig 2001 (RXTE, 1997) 


1.7±0.3 


7.4+0.2 


17.5+0.8 


4.7+0.2 


0.9(56) 


4.8 


2.5-30 


not reported 


Corbet & Peelel200ll(7?Xrg. 1997-1) 


1.71+0.03 


7.3+0.1 


17.3+0.6 


4.6+0.2 


0.82" 


3.1 


2-10 


not reported 


(RXTE, 1997-2) 


1.12+0.12 


5.3+0.2 


10.5+1.2 


2.7+0.7 


0.75° 


1.1 


2-10 


not reported 


Torrejon et sAXMZ\(RXTE. 1997) 


1.6+0.1 


7.6+0.4 


16.3+1.2 


4.5+0.4 


0.71(52) 


2.7 


2-10 


not reported 


(RXTE, 2001) 


1.6+0.1 


4.3+0.3 


20+2 


4.6+0.1 


1.27(49) 


1.3 


2-10 


29* 


(BeppoSAX, 1998) 


1.0+0.2 


7.8+0.5 


11+3 


1.1+0.3 


1.32(113) 


0.4 


2-10 


35+5* 


Masetti et al.l2004l(fleppoSAY. 1998) 


95+0" 


4 1+0.6 
^"'-0.5 


10.622 


0.8833 


1.1(219) 


0.4 


2-10 


35±5 C 


This paper (INTEGRAL, 2003, Rev. 67) 


1.8+0.7 


13+5 


22+6 


1.0 (fixed) 


1.2(154) 


15.9 


4-150 


32+5 


This paper (INTEGRAL, 2003, Rev. 87) 


1 7 +03 

1 ■ ' -0.4 


11+5 


29 + * 


1.0 (fixed) 


1.0(153) 


8.3 


4-150 


32+3 



" No information about the DOF is reported in the reference. * No significance reported. c At 2cr confidence level. 



bmc in xspec notation (Shrader & Titarchuk 1999), which pro- 
vided a ^^ ed of 1 .7 for 6 degrees of freedom (see top panel 
in Fig. |6j. We added a CSRF absorption feature, by using the 
xspec cyclabs model, at 32+2 keV to the bmc model, and ob- 
tained a slightly improved fit with of 1 .3 for 5 degrees of 
freedom. Finally, we also fitted the data by adding a Lorentzian 
profile in absorption to the bmc model. In this case the Xjied 
lowered to 1 . 1 for 5 degrees of freedom (see bottom panel in 
Fig-HJ. The center of the line was located at 31.5+0.5 keV and 



its FWHM was found to be 0.015+0.005 keV. The normaliza- 
tion of the line, 1.5+{JgXl0~ 3 count s~' keV -1 at a 68% confi- 
dence level (1<t), yields a significance of the line of ~2a J . We 



7 The significance of the CRSF detection could not be derived from 
the individual analysis of revolutions 67 and 87, due to the poor im- 
provement of the statistics obtained when adding the CRSF compo- 
nent in both cases. We can only state the significance of the detection 
after the nice improvement in the signal-to-noise ratio achieved by the 
mosaic. 
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Fig. 4. Spectra of 4U 2206+54 from JEM-X and ISGRI de- 
tectors. Top: The spectrum shown is for 2.2 ks of exposure 
time during revolution 67 (MJD 52761.36). Bottom: The spec- 
trum shown is for 6.6 ks of exposure time during revolution 87 
(MJD 52821.17). In both cases the spectral model shown in 
Table [2] is represented by the solid line, and the residuals to 
the model are displayed in the lower panels. The presence of 
an absorption feature around 32 keV, indicated by arrows, is 
suggested by the two datasets. 



note that the significance of this detection might be somehow 
influenced by changes in the continuum during the different 
observations. On the other hand, although the significance de- 
pends as well on the model chosen to fit the continuum, any 
model that properly fits the obtained continuum will reveal the 
presence of a ~2cr absorption around 32 keV, as can be seen 
from the data of Fig. [6] Thus we can conclude that the presence 
of a CRSF is strongly suggested by the data. This result, when 
combined with the previous claims after BeppoSAX and RXTE 
data (Torrejon et al. 2004 ; Masetti et al. 2004 1 gives evidence 
for the presence of this CSRF in 4U 2206+54. 

3.2. Radio 

No radio emission at 8.4 GHz was detected, with a 3<x upper 
limit of 0.042 mJy on 2003 May 12 (MJD 52771.4) and a 3cr 
upper limit of 0.066 mJy on 2003 May 20 (MJD 52779.7). We 



-4— 



4-U 



□ Crab 

• 4U 2206+54 



20 30 40 50 60 70 

Energy (keV) 

Fig. 5. Observed ISGRI spectra of the Crab and 4U 2206+54 
during revolution 67 (top panel), their ratio to the powerlaw 
model (middle panel), and the ratio of 4U 2206+54 to that of 
the Crab for the former ratios. The variations above 40 keV 
are not smooth, but random, and show large errors, indicat- 
ing that they reflect most likely extraction problems due to 
the fact that we are approaching the sensitivity limit of ISGRI 
(~0.04 count s _1 keV -1 for 2 ks exposures and 3cr level detec- 
tion at 45 keV). 

concatenated all the data and obtained a final 3<x upper limit of 
0.039 mJy. The resulting image is shown in Fig. [7] 

4. Discussion 

4.1. Conflicts in the neutron star scenario 

The aperiodic variability of the X-ray emission from 
4U 2206+54 favours the idea that the X-ray source is powered 
by wind-fed accretion onto a compact object. The X-ray lu- 
minosity of the source combined with its X-ray spectral shape 
likely excludes the possibility of a white dwarf (see, e.g., de 
Martino et al. 2004 for accreting white dwarf X-ray spectra). 
The compact object must be, therefore, a neutron star or black 
hole. 

There are two main difficulties for accepting 4U 2206+54 
as a typical wind-accreting neutron star. The first one is the 
lack of pulsations, as most other wind-fed systems are X-ray 
pulsars. In principle, this might result from a geometrical ef- 
fect: if the angle between the spin axis and the magnetic axis 
of the neutron star is close to zero or the whole system has a 
very low inclination angle, all the high-energy radiation seen 
could be coming from a permanently observed single pole of 
the neutron star. The system inclination is unlikely to be very 
small, as the projected v sin i for the optical companion is not 
particularly small (NR01), unless there is a very strong mis- 
alignment between the rotation axis of the optical star and the 
orbit. However, if the angles between the spin and magnetic 
axes of neutron stars are drawn from a random distribution, 
there is a non-negligible chance that for some systems they 
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Fig. 6. Top: Spectrum of 4U 2206+54 extracted from the mo- 
saicking of all images with significant detections from ISGRI. 
The solid line represents the bmc model spectrum (in xspec no- 
tation) fitted to the data. An absorption feature around 32 keV 
is seen in the residuals. Bottom: The same extracted spec- 
trum, but here the solid line represents the bmc model plus a 
Lorentzian feature in absorption fitted to the data, yielding a 
line position of 31.5+0.5 keV. 



will be aligned. Similar scenarios have been proposed to ex- 
plain the absence of pulsations from 4U 1700-37 (White et al. 
1 1 9 8 3 i> and also from the low-mass X-ray binary 4U 1700+24 
(Masetti et al. 2002 1, though in the latter case there is no con- 
clusive evidence that this system is sufficiently young to show 
pulsations. 

The second, stronger argument is the expected X-ray lumi- 
nosity. We can consider a canonical neutron star with 10 km 
radius and 1.4 M Q accreting from the fast wind of a low- 
luminosity 09III-V star in a 9.6 d orbit. Following the Bondi- 
Hoyle approximation, the accretion luminosity obtained, which 
is an upper limit to the X-ray luminosity, is ~lxl0 34 erg s _1 
(see Reig et al. 2003 for details about the method). In con- 
trast, the observed X-ray luminosity of 4U 2206+54 is in the 
range ~10 35 -10 36 erg s _1 , therefore comparable to those of 
HMXBs with OB supergiant donors (see Negueruela 2004|and 
references therein), which are believed to have mass-loss rates 
more than one order of magnitude higher than a 09 III-V star 
(Leitherer 1988 Howarth & Prinia [T989l . In addition, we note 
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Fig. 7. Image around 4U 2206+54, marked with a cross, ob- 
tained with the VLA at 8.4 GHz after concatenating data from 
2003 May 12 and 20. The image size is 4'x4'. Contours are -3, 
3, and 5 times the rms noise level of 0.013 mJy beam -1 . The 
ellipse in the bottom left corner represents the Full Width Half 
Maximum of the obtained synthesised beam of 10.6"x9.4" in 
PA=-83.4°. 



that our estimate for the semimajor axis of 4U 2206+54, of 55- 
60 R Q , is comparable to the highest values of semimajor axes in 
supergiant systems (see Kaper et al. 2004 1. Therefore, a close 
orbit cannot be invoked to solve the problem of the high X-ray 
luminosity. 

4.2. Excluding the black hole scenario 

If a black hole was present in the system, the photon index 
and luminosities of our INTEGRAL observations would indi- 
cate that the source is in a low/hard state (see McClintock & 
Remillard 2004 and references therein). Our radio observations 
took place on 2003 May 12 and 20, or during INTEGRAL revo- 
lutions 70 and 73, i.e., right between INTEGRAL observations 
at revolutions 67 (2003 May 01-04) and 87 (2003 Jun 30-Jul 
03). 

Gallo et al. (2003 1 found an empirical correlation between 
the soft X-ray flux (in the range 2-1 1 keV) and the centimetre 
radio emission (with observed flat spectrum in the range 4.9- 
15 GHz) for black hole binary systems in the low/hard state, of 
the form: S ra dio = (223 ± 156) x (5x) +0 ' 7 , where 5 ra dio is the 
radio flux density scaled to 1 kpc, Sx is the X-ray flux in Crab 
units scaled to 1 kpc, and the uncertainty in the multiplying 
factor is the non-linear lcr error of their fit. Therefore, by us- 
ing a measured X-ray flux we can compute the expected radio 
emission of a source in case it is a black hole. 

We obtained the flux from 4U 2206+54 in the 2- 
11 keV band from our JEM-X data, being 7.2 and 4.0xl0 10 



erg s cm for the 2003 May and June observations, re- 
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spectively. This flux was translated to Crab units by mea- 
suring the flux from the Crab in the 2-11 keV band, using 
JEM-X data from an INTEGRAL observation close in time to 
our 4U 2206+54 pointings. The Crab flux was found to be 
1.8 X 10~ 8 erg s^ 1 cm 2 , leading to fluxes of 40 and 22 mCrab 
for 4U 2206+54 during revolutions 67 and 87, respectively. 
From this, and using Nn = 1.0 xlO 22 atom cm (average of the 
values obtained by Torrejon et al. 2004 and Masetti et al. 2004 
from BeppoSAX data) we computed the unabsorbed corrected 
flux following equation (1) of Gallo et al. (2003 1, and then the 
resulting flux in Crab units scaled to 1 kpc distance (assum- 
ing a distance of 3 kpc to 4U 2206+54). The relation discussed 
above then predicts a radio flux density, already scaled back 
again to 3 kpc, of 12.6+8.8 mJy at the time of revolution 67 
and of 8.3+5.8 mJy for revolution 87 (where the errors come 
directly from the lcr uncertainties given in Gallo et al. 2003 
for the parameters of their fit). Thus, for revolution 67 the ex- 
pected radio emission would be in the range 3.8-21.4 mJy, and 
for revolution 87 in the range 2.5-14.2 mJy. 

We note that the lower expected radio flux density of 
2.5 mJy is already more than 60 times greater than the 
0.039 mJy 3<x upper limit found with our VLA observations. 
Obviously, it can be argued that our observations were not si- 
multaneous. During revolution 70 (which was in coincidence 
with the first radio observations) the flux found from ISGRI 
data is not significant enough and, unfortunately, the source 
was outside the FOV of JEM-X, which could have provided 
an X-ray flux suitable for this analysis. Nevertheless, we point 
out that the RXTE/ ASM count rate during our first radio ob- 
servation is very similar to that measured during revolution 
87, so it is reasonable to compare the obtained 2.5 mJy limit 
with our measured 0.042 mJy 3cr upper limit on that day, giv- 
ing again a difference of a factor ~60. It could also be ar- 
gued that the source could have experienced a transition to the 
high/soft state, that would naturally prevent the detection of ra- 
dio emission. However, in such a case the RXTE/ ASM count 
rates should have increased considerably, while during both ra- 
dio observations the count rates were similar or lower than dur- 
ing revolutions 67 and 87, when the photon indexes were typ- 
ical of low/hard states. In summary, if the correlation between 
X-ray emission and radio emission reflects indeed a general 
property of black hole systems, we conclude that there is not a 
black hole in 4U 2206+54. 

Moreover, Fender & Hendry ( 2000 ) show that all Galactic 
persistent black holes have detectable radio emission. As 
4U 2206+54 is a persistent source and does not show any de- 
tectable radio emission, it cannot host a black hole. Systems 
containing magnetised neutron stars (B > 10 11 G), on the other 
hand, do not show detectable radio emission. 

4.3. The cyclotron feature 

The indication of a cyclotron feature centred at 32 keV strongly 
suggests that there is a magnetic neutron star in 4U 2206+54, 
in good agreement with the lack of radio emission. INTEGRAL 
is the third mission reporting the likely detection of this absorp- 
tion feature (see Table|2j and, even if none of the detections can 



be considered statistically significant, the fact that it appears in 
three independent datasets cannot be ignored. 

If the line is indeed a CRSF we can compute the value of the 
magnetic field in the scattering region by means of the equation 
[fi/10 12 G] = [Soyci/ll -6 keV] (1 + z), being z the gravitational 
redshift at which we see the region. Considering that the line is 
produced at the surface of a canonical neutron star of 1 .4 M 
with a radius of 10 km, the gravitational redshift amounts to 
z=0.3 (see, e.g., Kreykenbohm et al. 2004), and from the po- 
sition of the line centre at 32 keV we obtain a magnetic field 
of 3.6 x 10 12 G. This value, in agreement with those found 
by Torrejon et al. (2004i and Masetti et al. (2004), is typical 
of magnetic neutron stars, and well within the range of 1.3- 
4.8 x 10 12 G obtained by Coburn et al. d2002> for a sample of 
ten X-ray pulsars displaying CRSFs (see their table 7). 

One is led to the conclusion that 4U 2206+54 is the first 
system known in which an accreting magnetic neutron star does 
not appear as an X-ray pulsar. In principle, the possibility of 
very slow pulsations cannot be discarded. The wind-accreting 
X-ray source 2S 01 14+650, with a Bl supergiant donor (Reig 
et al. I1996> . shows pulsations with a period of ~2.8 h (Finley 
et al. 11994b . The orbital period of the system is ~12 d, simi- 
lar to that of 4U 2206+54. Our timing analysis rules out the 
possibility of significant pulsations from 4U 2206+54 in the 
range 0.5-3.0 h. A modulation at a period of several hours is 
still possible, as existing datasets do not constrain strongly this 
period range. However, it seems more logical to conclude that 
geometrical effects are responsible for the lack of pulsations. 

5. Conclusions 

We present the first INTEGRAL GPS results on 4U 2206+54 to- 
gether with contemporaneous VLA observations of the source. 
A broad high-energy spectrum (4-150 keV), joining JEM-X 
and ISGRI data, has been extracted and fitted with spectral 
models similar to those used in the analysis of previously pub- 
lished data obtained with other satellites. The evidence for the 
presence of a cyclotron line is furthered, as INTEGRAL be- 
comes the third high-energy mission to observe a possible fea- 
ture at 32+3 keV. If the feature is indeed a cyclotron line, it 
indicates a magnetic field strength of 3.6 x 10 12 G, typical of 
magnetised neutron stars. 

Our VLA radio observations fail to detect the source at a 
very low level, indicating that any possible radio emission is 
at least 60 times weaker than what would be expected from 
a black hole system in the low/hard state. This lack of radio 
detection is again compatible with the presence of a magnetised 
neutron star. 4U 2206+54 appears to be the first known system 
containing an accreting neutron star that does not show up as 
an X-ray pulsar, most likely due to a simple geometrical effect. 

Longer high-energy exposures, such as an INTEGRAL long 
pointing, are needed to improve the S/N ratio to a level that will 
allow the confirmation of the presence of the cyclotron feature 
and the study of its eventual profile and energy changes with 
time and luminosity. Such observation would also allow the 
search for very long (several hours) or very weak pulsations. 
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